Measuring the Penetration Depth Anisotropy in MgB 2 using Small-Angle Neutron 

Scattering 
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Using small-angle neutron scattering we have measured the misalignment between an applied field 
of 4 kOe and the flux-line lattice in MgB2, as the field is rotated away from the c axis by an angle 
6. The measurements, performed at 4.9 K, showed the vortices canting towards the c axis for all 
field orientations. Using a two-band/two-gap model to calculate the magnetization we are able to 
fit our results yielding a penetration depth anisotropy, 7 a = 1.1 ± 0.1. 



PACS numbers: 74.25.Qt, 74.70.Ad, 61.12.Ex, 74.20.- 



The two-band/two-gap nature of superconductivity in 
magnesium diboride (MgB2) is well established^ In this 
material supercarriers reside on two distinctively differ- 
ent parts of the Fermi surface, with the larger energy 
gap (A CT ?s 7 mcV) originating from a bonding of p xy 
boron orbitals, and the smaller gap (A ff = 2.2 meV) 
arising from it bonding of the p z orbitals^ 3 ** Further- 
more the 7r band is isotropic, while the a band is nearly 
two-dimensional. 

The anisotropy of a typc-II superconductor is de- 
scribed either by 7 a = A c /A a b, where A is the magnetic 
penetration depth, or by j H = H c ^ ab / H c2 , c = £ab/€c, 
with £ being the coherence length. Traditionally these 
two anisotropies have been considered to be identical, 
but in materials with anisotropic gaps this is generally 
not the casei^i Magnesium diboride represents an ex- 
treme case where 7# ^ 7a. Theoretical work for this 
material indicates 7a (0) ~ 1.1 and 7h(0) w 6 at low 
temperatures, merging at J\(T C ) = jh(T c ) m 2.6 at the 
critical temperature T e i 5 i 7 i 8 

While there is experimental consensus on the temper- 
ature dependence of 7 g! 9 i 10 i 11 i 12 i 13 i 14 i 15 it has recently 
been suggested that superconductivity in MgB2 can be 
described by a single anisotropy factor, j(H) = jx(H) = 
7£(if), changing from j(H = 0) ~ 1 to j(H = H C 2) ~ 
6. 14,16 However, measurements of "j\{H) are still contra- 
dictory, with flux-line lattice imaging and torque magne- 
tomctry measurements yielding results ranging from an 
arguably field-independent value of 7a ~ 1.2^^^ to a 
strongly field-dependent 7a increasing with field and at- 
taining a value of ~ 3.5 at a modest field of 5 kOei 1 ^ " 
Measurements of the electronic specific heat in MgB2, 
which mainly reflects 75, yields an anisotropy which in- 
creases with the applied ficldi^i. 

In this paper we will demonstrate an use of small-angle 
neutron scattering (SANS) to determine the penetration 
depth anisotropy in MgB2, by measuring the misalign- 
ment between the applied field, 7i, and the direction of 
the vortices in the flux-line lattice (FLL) as the field is 



rotated between the c axis and the basal plane. Using 
a two-band/two-gap model we can fit the angular de- 
pendence of the misalignment. At 4.9 K and 4 kOe we 
find 7a = 1.1 ± 0.2, smaller than SANS measurements 
of the FLL anisotropy in single crystals^ and torque 
magnetometrj*i&, but in good agreement with SANS re- 
sults on powdersii and scanning tunneling spectroscopy 
results^ The misalignment corresponds to vortices cant- 
ing towards the crystalline c axis. 

The experiments were performed at the NG3 SANS 
instrument at the NIST Center for Neutron Research. 
The sample was a ~ 200 fig platclikc single crystal with 
the c axis parallel to the thin direction, grown using iso- 
topically enriched n B to reduce neutron absorption^. 
The n B concentration in the sample is estimated to be 
~ 95%, with ~ 5% 10 B due to the use of a natural 
boron BN crucible (20% 10 B, 80% n B) in the crystal 
growth. Incident neutrons with a wavelength A„ = 1 
nm and a wavelength spread AA„/A„ = 14% were used, 
and the FLL diffraction pattern was collected by a 650 
mm x 650 mm position sensitive detector with 5 mm 
resolution. Both the collimation and sample-to-detector 
distances were 8 m. All measurements were performed 
in a magnetic field of 4 kOe applied at different angles, 
6 = 0°— 80°, with respect to the crystalline c axis. For 
all angles the field was rotated around the vertical crys- 
talline [110] axis (see Fig. 1(a)), and the measurements 
were obtained following a field cooling procedure to a 
temperature of 4.9 K. The applied magnetic field is in 
the horizontal direction and approximately parallel to the 
direction of the incident neutron beam. 

Fig. 1 shows FLL diffraction patterns obtained for 
three different orientations of the applied field. Each im- 
age is a sum of the scattering from the FLL as the cry- 
omagnet is rotated around the vertical axis in order to 
satisfy the Bragg condition for the different reflections. In 
Fig. 1 (a) the field is applied parallel to both the sample c 
axis (9 = 0°) as well as the direction of the neutron beam. 
Consistent with our earlier reports a hexagonal FLL is 
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FIG. 1: (Color online)FLL diffraction patterns for MgB2 obtained at 4 kOe and 4.9 K, and with the applied field at angles 
9 — 0° (a), 45° (b), and 80° (c) with respect to the c axis. The orientation of the sample crystallographic axes is shown in 
panel (a). The reflections on the vertical axis in panel (c) are absent since they do not satisfy the Bragg condition. The data 
are smoothed by a 5 x 5 box average. 



observed, with the lattice plane normals (Bragg peaks) 
along the a axis^ Applying the field at an angle with 
respect to the c axis leads to a distortion of the hexagonal 
FLL as shown in Fig. 1(b), and eventually a reorienta- 
tion shown in Fig. 1(c). In contrast to H. \\ c where the 
reorientation progresses continuouly above an onset field 
H* ~ 5 kOe^S the transition in rotated fields is first or- 
der and occurs at 8 ~ 70° with an applied field of 4 kOe. 
The FLL reorientation is attributed to the multigap na- 
ture of the superconductivity in this material £L2& Pre- 
viously the FLL anisotropy has been used as a measure 
of 7^«i2i2P. However, recent measurements have indicated 
that in MgB2 in the proximity of the reorientation tran- 
sition this anisotropy may be influenced by vortex core 
effects^ similar to what is observed in members of the 
borocarbide superconductors^ In contrast to the FLL 
anisotropy, the direction of the vortices in the FLL does 
not depend on vortex-vortex interactions in the plane 
perpendicular to the flux lines. Hence measurements of 
the misalignment of the FLL with respect to the applied 
field provides a method for determining 7^, which is not 
affected by vortex core effects. 

The exact orientation of the applied magnetic field is 
obtained from the FLL rocking curve with H \\ c. In 
this case the vortices are parallel to the applied field. In 
Fig. 2(a) the intensities of the Bragg peaks lying on the 
horizontal axis are plotted as a function of the cryomag- 
net rotation angle, a. The calibrated zero (a = 0) is 
determined from the position of the midpoint between 
the two peaks, obtained by Gaussian fits to the rocking 
curves. As the field is rotated away from the c axis the 
vortices in the FLL cant away from the direction of the 
applied field as shown in Figs. 2(b) and 2(c). The rocking 
curve midpoint now directly reflects the misalignment of 
the FLL, 69, with respect to the direction of the applied 
field. Note that above the FLL reorientation transition 
four Bragg peaks shown in Fig. 1(c) are centered around 
the horizontal axis and used to determine 69. The mis- 
alignment corresponds to the vortices being rotated to- 
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FIG. 2: FLL rocking curves for MgB2 with the applied field 
at 8 — 0° (a), 45° (b), and 80° (c) with respect to the c axis. 
The intensities at each angular setting, a, is obtained by sum- 
ming the detector counts at position of the Bragg reflections. 
The curves are Gaussian fits to the data. No background 
subtraction is performed. 
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FIG. 3: FLL misalignment (58) as a function of the applied 
field rotation (8). The error bars take into account both the 
statistical and systematic errors of the measurements. The 
solid line is a fit to the data described in the text. 



xz-plane forming an angle with z. The linear relation 
between the applied field H, the induction B, and the 
magnetization M reads' 



Bi 



Hi + 4irMi(l — n j) , i = x, z , 



(1) 



where n z = 1 — 2n and n x = n are the demagnetization 
coefficients. For a coin-shaped sample n<l. The angle, 
(9b = 9 — 69, between the induction and the c axis is given 
by 

, a B x H x + 4wM x (l - n x ) 

tant^R = — = ; f . (2) 

B z H z + 47rM z (l-n z ) 1 ; 

For Tt H c i where M <C H, as is the case in the present 
experiment, the misalignment 59 = 9 — 9g is thus given 
by: 



wards the crystalline c axis. Only a slight increase in the 
average rocking curve width from 0.45° ± 0.05° FWHM 
(0 = 0°) to 0.55° ± 0.05° FWMH (0 = 80°) is observed 
as the field is rotated away from the c axis, indicating 
that no significant disordering of the FLL takes place. 
The experimental resolution is estimated to be ~ 0.25° 
FWHM. The absolute scattered intensity from the FLL 
decreases due to the increased absorption, as the effective 
sample thickness becomes larger. 

The results of the SANS measurements are summa- 
rized in Fig. 3. This shows a symmetrical 69 vs. 9 
dependence, with a maximum 59 ~ 1° close to 9 = 45°, 
and decreasing towards zero for 9 = 90° as expected. 
Measurements at angles above 80° were not possible due 
to absorption and increasing background scattering from 
sample defects. 

While the torque in a single-band superconductor is 
proportional to the misalignment, 59, this is not the case 
is materials with 7a ^ jh such as MgB2 . In the following 
it is assumed that the sample can be approximated by a 
oblate spheroid (a coin) with the c axis along z, and 
choosing (x, y) so that the applied field TL lies in the 



2- 



= sin 20 



M z M x , 
rtz rL x 



(3) 



For materials with 7^ 7^ 'Jh the free energy density is 
given by£ 



B 2 

F = F — — — 



n ^ c7 it , (4) 



8t 32^ A 2 7 y 3 Be x (l + (3)^ 
where A = (A^Ac) 1 / 3 is the average penetration depth, 



7a h cos b , and 




7a £h 
1h£\ 



(5) 



The constant r\ lumps together the London uncertainty 
factor ?s 1.4 with yielding 7/ f» 2.8. 

Evaluation of the magnetization, M = —dF/dB, and 
the misalignment angle, 50, is now straightforward. Dif- 
ferentiating with respect to (B X ,B Z ), and noting that 
Be x = y/Bl+Hl 
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where Mo = <f>o / S2tt 2 \ 2 . Since M -C ~H, the angle 50 approximation. Substituting Eqs. JfjJ and into Eq. 
is small and Mi's in Eq. © can be taken in the lowest © one can replace (B,0b) with (Ti.,0) obtaining: 
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For MgB2, 7h > 7a at all temperatures, and hence 
89 > for all values of 9 and T provided that n is small. 
Consequently, for a coin-shaped sample, the vortices are 
expected to rotate towards the c axis, consistent with the 
experimental results. 

The solid line in fig. 3 is the best fit of Eq. (jSJ to the 
SANS data, obtained by varying 7a and A, while setting 
"/H = 6, r) = 2.8 and using an approximate demagnetiza- 
tion factor for the sample, n = 0.01. This yields a value 
7a = 1.1±0.2 and A = 155 nm. The results provides addi- 
tional evidence for a low value of 7a even at this relatively 
high field. Due to the uncertainty on the parameter r\ as 
well as the demagnetization factor, n, a reliable determi- 
nation of A and 7# is difficult, especially since changing 
these two parameters both affects only the amplitude of 
89. Nonetheless, the fitted value of A is within the range 



of literature values summarized in Ref. ||| and in fair 
agreement with A = 136 nm from Ref. |l6j |. However, we 
can not exclude a field dependent 7^ (H ) merging with jh 
at H c2 , as suggested in Rcfs. [14, 16, 21]. On the other 
hand the profile of 89 is determined solely by 7a, allowing 
a reliable fit of this parameter. Increasing 7a causes the 
<5(9-profile to become more asymmetric as shown in Fig. 
4(a). Although it would be possible to fit our data with 
7# = 7a = 1-1 and A = 135 nm, this is an unreasonably 
low value of jh for a field of 4 kOe, and consequently the 
results presented here do no t supp ort a single anisotropy 
factor as proposed in Refs. |l4ll6l |. 



Finally, it is of interest to compare 89 to the torque 
density, r y = M z Ji x — M x 7i z . Inserting the Mi's from 
Eqs. JfjJ) and one obtains^ 
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Comparing this to Eq. J3J) it is clear that in general 
89 is not proportional to r y as only the former is shape 
dependent. In particular Eq. 10 suggests that at low 
temperature the torque is negative as shown in Fig. 
4(b), in contrast to what has been observed for layered 
superconductors^^ Near T c one recovers 89 tx r y as in 
materials with 7a = 7h . Given the sucsess of the model 
presented here in describing the misalignment, it would 
be of interest to test further the relationship between 89 
and the torque. Finally, one notes that for a spherical 
sample with n = 1/3 one finds 89 oc t v for all values of 
7a and j H - 

In summary we have shown that SANS measurements 
of the misaligment angle between the applied field and 
the FLL can provide information about the penetration 
depth anisotropy in MgB2. At 4.9 K and 4 kOe we find 



a relatively low value of 7a = 1.1 ± 0.2. More studies are 
needed to extend the measurements of 7a to both higher 
temperatures and fields, in order to investigate the entire 
mixed phase of MgB2 . 
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